Incommensurability and spin excitations of diagonal stripes in cuprates 
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Based on the time-dependent Gutzwiller approximation we study the possibiUty that the diag- 
onal incommensurate spin scattering in strongly underdoped lanthanum cuprates originates from 
antiferromagnetic domain walls (stripes). Calculation of the dynamic spin response for stripes in 
the diagonal phase yields the characteristic hour glass dispersion with the crossing of low energy 
Goldstone and high-energy branches at a characteristic energy Ecross at the antiferromagnetic wave- 
vector Qaf- The high energy part is close to the parent antiferromagnet. Our results suggest that 
inelastic neutron scattering experiments on strongly underdoped lanthanum cuprates can be under- 
stood as due to a mixture of bond centered and site centered stripe configurations with substantial 
disorder. 

PACS numbers: 74.25.Ha, 71.28.-|-d, 71.45.Lr 



The low energy spin response of many high-Tc su- 
perconductors is characterized by magnetic fluctuations 
which are incommensurate with respect to the antifer- 
romagnetic (AF) order (for a review see [J]). Since 
in codoped (with Eu, Nd or Ba) lanthanum cuprates 
(LSCO) an accompanying charge-density wave has been 
detected (e.g. [^) this kind of scattering is usually at- 
tributed to the formation of stripe textures where the 
doped holes segregate into quasi one-dimensional 'rivers 
of charge' which simultaneously constitute domain walls 
for the AF order parameter. 

The similarity of low-energy inelastic neutron scatter- 
ing (INS) data between codoped and non-codoped mate- 
rials supports the picture of a fluctuating and (or) glassy 
stripe phase in LSCO. From these experiments it turns 
out that above doping Uh = 0.055 the stripes are ori- 
ented along the Cu-0 bond direction ('vertical stripes') 
and the incommensurability 6 (deflned as the deviation 
of the magnetic peak from Qaf) linearly increases up to 
nil ~ 1/8 [3]. Above rih « 1/8, S stays essentially con- 
stant but the intensity of the low energy spin fluctuations 
decreases and vanishes at the same concentration where 
superconductivity disappears in the overdoped regime 
thus suggesting an intimate relation between both phe- 
nomena [4|- 

Upon lowering the doping below rihj^ 0.055 the incom- 
mensurability d rotates by 45° 0, @, 0, Hi to the diagonal 
direction and the orthorhombic lattice distortion allows 
one to conclude that the elastic diagonal magnetic scat- 
tering is one-dimensional with the associated modulation 
along the orthorombic 5*-axis. When S is measured in 
units of reciprocal tetragonal lattice units in both the ver- 
tical and diagonal phase it turns out that the magnitude 
of the incommensurability numerically coincides across 
the rotation leading to a linear relation S — x. Upon 
approaching the border of the AF phase at — 0.02 
the incommensurability approaches e — x ^ where e is 
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FIG. 1: (Color online) Charge density and magnetization for 
site centered and bond centered diagonal stripes. The den- 
sities are plotted along an horizontal scan, indicated by the 
dashed arrow. 



measured in units of reciprocal orthorombic lattice units, 
thus e — \/2S. 

In this paper we present computations of the incom- 
mensurability and the dynamic response of diagonal 
stripes based on the Gutzwiller approximation and its dy- 
namical extension and compare with our previous compu- 
tations for vertical stripes and recent INS data [Ifl] . Our 
results suggest that in the spin glass phase incommen- 
surate scattering is due to a disordered mixture of bond 
centered and site centered diagonal stripe configurations. 
Calculations are based on the one-band Hubbard model 
(on-site repulsion U) with hopping restricted to nearest 
(~ t) and next nearest (~ t') neighbors. We first treat 
the model within an unrestricted Gutzwiller approxima- 
tion (GA) which leads to the stripe textures On top 
of this we obtain the spin excitations within the time 
dependent Gutzwiller approximation [13, [Sj- Parame- 
ters are U/t = 7.5, t' /t = -0.2 and t = SAOmeV^. 
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In previous papers we have shown that our approach 
can reproduce both, the magnon excitations of undoped 
LC0[3| as revealed by neutron scattering [l^ and the 
doping dependent incommensurability in the vertical 
stripe phase H, 17|. In addition, the formalism leads to 
quantitative agreement with the dispersion of spin excita- 
tions in and can reproduce the doping 
dependence of the optical conduct ivitv[20l|. 

Fig. [T] shows the charge and magnetization profile for 
the diagonal stripe structures we have obtained as stable 
saddle-point solutions within the GA. The first is a site 
centered diagonal (DSC) stripe where the doped holes 
are confined to nonmagnetic sites which simultaneously 
constitute the antiphase boundary for the AF order pa- 
rameter. 

The second stable texture (Fig. [TJ right panel) can be 
considered as the diagonal counterpart to vertical bond- 
centered (VBC) stripes. However, due to the diagonal 
orientation the stripe acquires a net ferromagnetic mo- 
ment. Moreover, the magnetization points in the same 
direction (i.e. the whole layer becomes ferromagnetic) 
when the BC stripes are separated by an odd number of 
lattice constants in x- (or y-) direction. They can be con- 
sidered as the smallest 'staircase' variant as discussed by 
Granath 2l|. We have found that more extended 'stair- 



case' structures are higher in energy and do not yield 
stable magnetic excitations. 

The approximate linear relation between incommen- 
surability and filling at low doping [3, d, H, 0, S| implies 
that isolated stripes are self-bounded linear aggregations 
of holes with a well defined number of holes per Cu along 
the stripe given by i^opt = nh/e- Indeed as for vertical 



stripes [16|, |17[ we find an optimum filling, however, an 



important difference here is that Vopt depends substan- 
tially on the texture. DSC are insulating (the Fermi level 
falls in a gap) with i^opt ~ 1, thus e « n^. Diagonal BC 
(DBC) stripes are metallic and have an optimum filling 
at Vopt ~ 0.75 implying 5 ^ Uh [23 |. 

While DBC stripes are practically (accidentally) de- 
generate in energy with low doping VBC, the energy of 
the DSC texture is « 0.02t per hole above [2§|. We be- 
lieve that these small energy differences are not signifi- 
cant given the simplicity of the model. A precise deter- 
mination of the relative stability of the different phases 
would require at least multiorbital effects, inclusion of 
both long-range Coulomb interactions and coupling of 
the holes to the tilts of the CUO4 octahedra. The latter 
have been shown to play a major role in the stabilization 
of vertical vs. diagonal stripes [is^. On the other hand 
we believe that results within one phase are much more 
reliable. For example we have checked that the different 
optimum v for DSC and BSC stripes for a 3-band hamil- 
tonian do not differ from those found for the one-band 
model so that we consider this as a robust feature of our 
calculation. 

In the present case the more important extrinsic effect 
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FIG. 2: (Color online) Doping dependence of the magnetic 
incommensurability e in the diagonal phase as deduced from 
the referenced elastic neutron scattering experiments. Solid 
and dashed line are the predicted dependence for DSC and 
DBC stripes, respectively. 



is disorder. For vertical stripes the correlation length can 
reach IbOatetra or around 20 times the magnetic stripe 
periodicity while for diagonal stripes the correlation is of 
the same order or even smaller than the periodicity 6] . 
This should be kept in mind while comparing with our 
results which correspond to perfectly ordered stripes ar- 
rays. Fig. [2] compares the doping dependence of the mag- 
netic incommensurability e in the diagonalphase as seen 
in elastic neutron scattering experiments 0, 0, S [l^l to- 
gether with the predicted dependence for DSC and DBC 
stripes. For non-codoped samples the experimental in- 
commensurability is in between the theoretical ones for 
DBC and DSC, with a stronger tendency for DBC close 
to the metal-insulator transition which shifts to DSC at 
low doping. An interesting effect occurs upon substitu- 
tion of Cu^+ by nonmagnetic Zn^+ which has a filled id}'^ 
shell 2J| . This should favor DSC stripes which are forced 
to have zero magnetic moments in the core and indeed 
the incommensurability shift to the computed DSC line 
(c.f. Fig. (2) is in accord with our expectation. All this 
supggests that the system is made of a disordered mix- 
ture of DSC and DBC configurations with a gradual shift 
from the former to the latter as doping is increased. 

In the following we discuss the spin excitations for 
doping Uh = 0.05. For DSC stripes we have e = 0.05 
(i.e. charged stripes separated by 20a along the x-, or 
y-direction and for DBC stripes with e = 0.071 (i.e. sep- 
arated by 14a along x-, or y-direction). Fig. O displays 
/(ijj,q) = u}S{uj,c\) with S{u!,c\) the transverse dynam- 
ical structure factor relevant for magnetic neutron scat- 
tering experiments. The uj factor makes visible the high 
energy features which would appear much faint in the 
experimentally accessible S{lo,c\). We show the disper- 
sions of magnetic excitations perpendicular (left panel) 
and parallel (right panel) to the stripe direction. 
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FIG. 3: (Color online) Magnetic excitations perpendicular 
and parallel to DSC (upper panel) and DBC (lower panel) 
stripes (doping Uh ~ 0.05) ad reveled by a contour plot of 
/(a;,q) = tjS'(aj,q). The momentum scale is the distance 
from Qaf measured in units of 2ti / aortho- 



Parallel to the stripe the dispersion is dominated by 
modes very similar to the ones of the undoped AF as 
could be expected from the real space structure. Due to 
the large unit cell and the small extension of the reduced 
Brillouin zone in the perpendicular direction excitations 
consist of a large number of folded bands. These are 
generally silent when plotted in the extended Brillouin 
zone but acquire a large spectral weight in proximity to 
the modes of the parent antiferromagnet (shown with a 
dashed line) giving rise to the horizontal segments in the 
right panel. This is analogous to the effect reported in 
photoemission where spectral weight is maximized when 
intersecting the free electron dispersion relation 2]J . This 
clarifies the evolution of the spectrum at very low dop- 
ing: As doping decreases the dynamical structure factor 
gradually tends to the one of the antiferromagnet by in- 
creasing the number of bands and modulating the spec- 
tral weight. 

At low energies the spectrum is determined by the 
Goldstone excitations belonging to an inwards (i.e. in the 
direction of [tt, tt]) and an outwards dispersing branch. 
For DSC stripes upon increasing energy the inwards dis- 
persing branch dominates in weight and reaches the AF 
wave-vector Qaf at Ecross ~ 35meV where it connects 
to the dispersion along the stripe. In contrast to the ver- 
tical excitations [ij] the saddle-point structure is very 
small and barely resolved due to finite size effects. For 
DBC stripes a small gap appears between the Goldstone 
mode and the saddle point in strong resemblance with 
linear spin- wave theory (LSWT) results by Carlson et al. 
[2^. There it was pointed out, that for even spaced di- 
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FIG. 4: (Color online) Constant frequency scans of Imxq'jJ 
for DSC stripes (rih ~ 0.05). Data have been convoluted 
by a gaussian in frequency and momentum space which half- 
width corresponds to the resolution of the calculation, i.e. 
Slo — O.Olt and Sq — 2n/ 801 /atetra- The center of each panel 
corresponds to Qaf and wave-vectors are denoted in units of 

2Tv/atetra- 



agonal stripes (as in the present case) Qaf becomes a 
reciprocal lattice vector with a downturn of the acous- 
tic branch. Interestingly for DSC stripes in our itinerant 
approach the analogous gap is not resolved resulting in 
practice in a continuous dispersion of the Goldstone mode 
up to the saddle with energy Ecross- We expect that in 
both cases the gap will be washed out by disorder and 
damping effects. 

In Fig. m we show cuts of the dispersion for the DSC 
stripes. At low frequencies the excitations at the two in- 
commensurate wave-vectors carry the dominant weight 
which start to form spin-wave cones upon increasing en- 
ergy. The cones merge at Ecross = 0.1 2t where it is ap- 
parent that the weight of the cones is strongly anisotropic 
and enhanced close to Qaf- Similar to our previous find- 
ings for vertical stripes 14, 1^ the response well above 
Ecross becomes two-dimensional forming a ring shaped 
pattern around Qaf- 

Fig.Oreports the doping dependence of Ecross for both 
SC (main panel) and BC (lower right panel) textures in 
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FIG. 5: (Color online) Doping dependence of Ecross for DSC 
and DBC stripes. Stars refer to results for vertical stripes 
and squares to diagonal ones. Circles with error bars are 
experimental data from Refs. [lol. [2^. [27l. [2^. Dashed lines are 
the extrapolation of Ecross to = at = 0. 



the vertical and diagonal phase compared to INS data 
lol . 2^ . 27 . 2^. Since for the undoped system one recov- 
ers the spin excitations of the AF, Ecross extrapolates 
down to zero frequency at n/j = (dotted line) for both 
DSC (main panel) and DBC (lower right panel) stripes. 
Due to the large supercells involved we could not access 
the doping regime Uh < 0.04. 

Below X = 0.08 our computations overestimate Ecross- 
One should keep in mind that we are working with a 
fixed set of parameters determined in Ref. 1 14 - For exam- 
ple reduction of t improves considerable the agreement 
but spoils the fitting of the undoped dispersion relation. 
Due to the complexity of the computations we have not 
performed intense optimization of the parameter set, so 
it may be that for different choices of U and t' and t 
we could also achieve a better fit of the overall doping 
dependence. However, we believe the main reason for 
the disagreement is the strong disorder character of the 
stripes in this doping range. Indeed quite generally one 
expects a softening of excitations as the magnetic config- 
uration becomes more disordered. This becomes evident 
in the Ising limit comparing the dispersion relation of 
an Ising ferromagnet (or AF) with that of an Ising spin 
glass in which the exchange constant J is the same but 
the sign is random. In the ordered case the spin flip exci- 
tation energy is 2zJ with z the coordination, whereas in 
the glassy case it is smaller in average due to frustrated 
configurations. We have checked for an array of stripes 
on a finite cluster that the inclusion of random local im- 
purity potentials leads to a softening of Ec oss, whereas 



a variation of the stripe spacing induces a broadening of 
the corresponding excitation. 

To conclude, we have shown that the doping depen- 
dence of the incommensurability of diagonal magnetic 



scattering is compatible with a mixture of SCD and BCD 
stripes. Assuming an ordered array of stripes we find 
excitations that are qualitatively in accord with experi- 
ment but too high in energy. We attribute the dif- 
ference to the much more disordered nature of diagonal 
incommensurate scattering. Our intermediate coupling 
approach predicts a quasi two-dimensional high energy 
response in the diagonal phase of LSCO which displays 
a one-dimensional static response along the orthorombic 
6*-axis. Our computations also clarify how the disper- 
sion relation converges to the one of the AF as doping is 
decreased. 
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